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Abstract 



Boundary value problems for integrable nonlinear partial differen- 
tial equations are considered from the symmetry point of view. Fam- 
ilies of boundary conditions compatible with the Harry-Dym, KdV 
and MKdV equations and the Volterra chain are discussed . We also 
discuss the uniqueness of some of these boundary conditions. 



1 Introduction 

In our previous paper [1] we have briefly discussed a method to construct 
boundary value problems of the form 

Ut = f(u,Ux,U2, — ,M n ), (1) 

p(u,u 1 ,u 2 ,...,u k )\ x=Q = 0, (2) 

completely compatible with the integrability property of Eq.(l). Here u = 
u(x,t) , Ui = |^ and / is a scalar (or vector) field. The aim of the present 
paper is to expound detailly our scheme and also extend it to the integrable 
differential-difference equations. 
Let the equation 

u T = g(u,ui,...,u m ), (3) 

for a fixed value of m, be a symmetry of the equation ([!]). Let us intro- 
duce some new set of dynamical variables, consisting of the variable v = 
(u, Mi, m 2 , ...u n -i), and its ^-derivatives v t , v u , ■■■ ■ One can express the 
higher x-derivatives of m, i.e., Mj for i > n and their ^-derivatives, by the 
utility of the equation ([]]), in terms of the dynamical variable v and their 
t-derivatives. Here n is the order of the equation (]]]). In these terms the 
symmetry (^|) may be written as 

v T = G(v,vt,v t t,...vtt...t). (4) 

We call the boundary value problem, Eqs.([T|) and (Q), as compatible with 
symmetry (Rf) if the constraint p(v) = (or constraints p a (v) = where 



a = 1,2,...N and N is the number of constraints) is consistent with the 
r-evolution 

-^ = 0, (mod p = 0) (5) 

or 

Eq.fl5|), by virtue of the equations in (f|), must be automatically satisfied. In 

fact (|^) means that the constraint p = defines an invariant surface in the 

manifold with local coordinates v. This definition of consistency of boundary 

value problem with symmetry is closer to the one introduced in [3], but not 

identical. For instance, let us examine whether the boundary value problem 

Ut = u xx \ u x = cu, x = 0, is compatible with the symmetry u T = u xxx . To 

this end one has to check if the equation w = cu defines an invariant surface 

for the system of equations u T = w t , w T = u u (here w is U\). Evidently 

the answer is negative. To check the validity of compatibility condition in 

the sense of [3] one has to compare two sets of equations u 2n +i = cu2 n and 

«3 n+ i = cu-in , n > 0. These equations are obtained by differentiation of the 

constraint equation u x = cu with respect to t and r variables respectively. In 

this sense the boundary condition is compatible with the symmetry because 

the sets of equations don't contradict each other. 

We call the boundary condition fl2|) is compatible with the equation if it 
is compatible at least with one of its higher order symmetries. 

Our main observation is that if the boundary condition is compatible with 
one higher symmetry then it is compatible with infinite number of symme- 
tries. We define a set S with infinite number of elements where its elements 
are symmetries of the Eq.(|l|). S may or may not contain the whole symme- 
tries of (|ip. For instance , S contains the even numbered time independent 
symmetries for the Burgers equation. 



We note that all the known boundary conditions of the form (|2]) consistent 
with the inverse scattering method are indeed compatible with the infinite 
series of generalized symmetries. On the other hand, stationary solutions 
of the symmetries compatible with @ allow one to construct an infinite di- 
mensional set of "exact" (finite gap) solutions of the corresponding boundary 
value problem (p]) and (0). However , in this work we do not discuss analyt- 
ical aspects of this problem. We note also that , in this paper we shall deal 
with boundary conditions of the form given in @. An effective investigation 
of boundary conditions involving an explicit ^-dependence is essentialy more 
complicated. Such a problem has been studied , for instance , in JHJ . 

The article is orginized as follows. In Section 2 we introduce the necessary 
notations and give the Burgers equation as an illustrative example. We prove 
that if a boundary condition is compatible at least with one higher order 
symmetry then it is compatible with every one of even order. In Section. 3 
we consider the nonlinear Schrodinger, Harry-Dym, Korteweg de Vries and 
modified KdV equations. Using the symmetry approach we find a boundary 
conditon compatible with the symmetry algebra of the Harry-Dym equation 
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u x = cu , x = (6) 

u xx = c 2 u/2 , x = 0, 
where c is an arbitrary real constant. Actually one has here two constraints. 
Although we are taking the boundary conditions at x = , one can shift this 
point to an arbitrary point x = Xq without loosing any generality. We con- 
jecture that the boundary value problem given in (U) is compatible with the 
Hamiltonian integrability and solvable by the inverse scattering technique. 
In addition we conjecture that (using the idea in [2]) one can prove that on 



the finite interval x\ < x < X2 the Harry-Dym equation with the boundary 
conditions u x = c$ u, u xx = Cq u/2 for x = x\ and u x = c\ u, u xx = c\ u/2 for 
x = X2 is a completely integrable Hamiltonian system. 

Section 4 is devoted to the differential-difference equations. In the last 
section we propose further generalization of the compatibility and discuss on 
some open questions. 

2 Boundary Conditions Compatible with Sym- 
metries 

In the sequel we suppose that eq. (|l[) admits a recursion operator of the form 

(see [4-6]) 

R = J2^D i + J2 a -^D- 1 a^ l , n>0 , h>0 (7) 

i=0 i=Q 

where ctj, a_i : i, a_2,; are functions of the dynamical variables, D is the 
total derivative with respect to x. Recursion operators when applied to a 
symmetry produce new symmetries. Passing to the new dynamical variables 
v,Vt,Vu, ■■■ one can obtain, from (|7]), the recursion operator of the system of 
equations (|j) (we don't prove that every recursion operator may be rewritten 
in the matrix form, but we will give below the matrix forms of the recursion 
operators for the Burgers, KdV, MKdV and Harry-Dym equations) 

M K 

R = ^a i (9 i ) i + ^a_ M (9 t - 1 )a_ 2il , M > 0, K > 0, (8) 

i=0 i=0 

where a; depends on v and on a finite number of its t-derivatives , dt is the 
operator of the total derivative with respect to t. If (P is a scalar equation, 



R is a scalar operator, then R is an n x n matrix valued operator. Our 
further considerations are based on the following proposition, which really 
affirms that if an equation admits an invariant surface, then a kind of its 
higher symmetries admits also the same invariant surface. 

Proposition 2.1. Let the equation (|]) is of the form v T = H(R)v t , 
where R is the recursion operator (8) and if is a polynomial function with 
scalar constant coefficients. If this equation is consistent with the constraint 
p(v ) = 0, where rank of p equals n — 1 (here n is the dimension of the vector 
v) then every equation of the form v T = L(H(R))v t , where L is arbitrary 
chosen scalar polynomial with constant coefficients, is also compatible with 
this constraint. 

Proof: Introduce new variables w = (w 1 , w 2 , ...w n ) in the following way: 
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W = p , W 



P 



.w 



n-1 



P 



n-1 



and w n = p n is a function of v, here p l 



is a coordinate of the vector p. Then one obtains the equation w T = Pwt 
from (5), where P = A~ 1 H(R)A and A = dv/dw is the Jacobi matrix of the 
mapping w — > v. Notice that under this change of variables the constraint 
p{v) = turns into the equation w l = for i — 1,2, ...n — 1. Imposing this 
constraint reduces the equation w T = Pwt to the form 

/ \ / P u ... P ln \ / \ 

Pra-1,1 ••• Pn-l,n 

\w n T ) V P»,l .- Pn,n J \<y 

Let us show that elements of the last column of the matrix P are equal zero 
except maybe P„ in : Pj, n = 0forl<i<n — 1. Really, letting P^ n ^ 
for some j < n — 1 the equation P^ n w™ = gives a connection between 
variables w n , w", ... which are supposed to be independent. Since the set 



of such block triangular matrices constitutes a subalgebra in the algebra of 
the all squared matrices hence one can easily conclude that the operator 
L(P)(modw l — 0, % < n — 1) is also block triangular so the equation w T = 
L(P)wt is consistent with the constraint w l = 0, % < n — 1). It completes the 
proof of the Proposition 2.1. 

Proving later a kind of uniquness theorem (see below the Proposition 2.2) 
we will use the following statement 

Proposition 2.1': Suppose that the constraint p(v) = of the rank equal 
n — 1 is compatible with the equation (5) having the form v T = H(FL n °)v t , 
where H = H(z) is a scalar polynomial function of z with constant coef- 
ficients. Assume that no > 1 is an integer and the leading term b^ (the 
coefficient before the highest derivative) in the expression R n ° = b^{d t ) N + 
^v-i(<9t) iV_1 + ••• is a scalar matrix, i.e. it is proportional the unit matrix. 
Then the constraint is consistent with the equation 

v T = R no v t (9) 

Proof: In terms of the variable w we have introduced proving the previ- 
ous proposition the equation v T = H(R n °)vt takes the form w T = H(Ri°)wt- 
Owing the fact that the point transformation preserves the commutativity 
property of flows , the operator Ri = AHA~ l is the recursion operator in the 
new variables. Again, just in the previous proposition one has that the oper- 
ator P = H(Ri°) under the substitution p = (or really, w l = 0, i < n — 1) 
is a lower block triangular matrix valued one. Our aim now is to prove that 
the operator Q = R"° (modw 1 = 0, % < n— 1) is also block triangular. Setting 
H(Q) = a n Q n + a n -\Q n ~ l + ... + a® and representing Q as formal series 



X)fcl_oo c kdf using the famous Campbell-Hausdorf formula one obtains that 

H(Q) = a n (c n M (d t ) nM + nc^ l c M ^{d t ) nM - 1 + ...) + ... + a 

One has that H(Q) belongs the set M_ consisting of the all lower block 
triangular matrices. By looking at the coefficients of different power of the 
operator dt one can show that the matrices q, % = M — 1, M — 2, ... satisfy 
the equations 

c^c. + T.eM., 

where Tj are polynomials with scalar coefficients on variables c i+ i, c i+2 , ■■■Cm 
and their derivatives. So, because of assumptions cm = ^m £ Mq and 
detbM 7^ it is easy to prove by induction that q G M for all i < M. 

For illustration let us give the Burgers equation as an example. The 
Burgers equation 

u t = u xx + 2uu x (10) 

which possesses the recursion operator of the form 

R = D + u + u x D- 1 (11) 

(see, for instance, [7]). The simplest symmetry of this equation is u T = u x . 
In terms of the new dynamical variables this symmetry equation takes the 
form 

(12) 
Ui iT = u t — 1UU\ 

This equation does not admit any invariant surface of the form p(u, Ui) = 0. 
Really, differentiating this constraint with respect to r one obtains 
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^« 1 + ^-(tt t -2«« 1 ) = (13) 

OU OU\ 

Because of independence of the variables u t and u\ we have 

£- = £ = <> (14) 

OU\ OU 

which leads to a trivial solution p = constant. As a conclusion we don't have 
any invariant surface (curve) in (u, Ui) - plane. Similarly the third order 
symmetry u T = u 3 + 3uu2 + 3u 2 + 3u 2 U\ rewritten in the new variables (u, U\) 
gives the following system of two equations 

u T = u l:t + uu t + (u 2 + u 1 )u 1 , 

(15) 

Ml,r —U tt — U Ui >t + (u 2 + Ui)u t -2u Mi (u 2 + Mi) 

This system also does not admit any invariant surface of the form 
p(u, Mi) = 0. It may be easily proved that the same is true for every symmetry 
of the odd order , i.e ., u T — U2 m +i + h(u2 m , ...,u). Because the correspondent 
system of equations has different orders in the highest t-derivatives 

Mr = ( 9 t m Ml + ... ? u ltT = d™ +l u + ... (16) 

Unlike the symmetries of odd order, for the symmetries of even order the 
correspondent system of equations has the same orders in the highest it- 
derivatives. This fact leads us to show that the symmetries of even order 
admit an invariant surface p(u,Ui) = 0, depending upon two arbitrary pa- 
rameters. 

Proposition 2.2: If the boundary condition p(u,Ui)\ x= o = is compat- 
ible with a higher symmetry of the Burgers equation, then it is of the form 
(see [3]) c(u\ + u 2 ) + C\ u + c 2 = and is compatible with every symmetry of 



the form u T = P(R 2 ) Ut where P denotes polynomials with scalar constant 
coefficients. 

Proof:The Frechet derivative of (|10D gives the symmetry equation of the 
Burgers equation 

d t a= (D 2 + 2uD + 2w)a (17) 

where w stands for u\. As the operators acting on symmetries we may take 

IT 1 = 0- 1 (D + 2u) (18) 

in the recursion operator flTTp. Consequently the recursion formula u n+1 = 
Ru Ti becomes 

u n +1 = (u + 2 w d^ 1 u) u Ti + (1 + w <9 t _1 ) w Ti (19) 

Differentiating it with respect to x and replacing w x = u<i = u t — 2 u w one 
obtains 



w n+1 = [d t + 2{u t -2u w) djT 1 u] u n + \-u + (u t -2uw) d^ 1 ] w n (20) 
for i = 1, 2, .... Thus the matrix form of the recursion operator R is given by 

R= * , ' _J ( 21 ) 

\ d t + 2 (u t — 2 u w) d t u —u + (u t — 2 u w) d t ) 

It is well known that every higher order local polynomial symmetry may be 
represented as a polynomial operator Pq(R) applied to the simplest classical 
symmtery u T = u x . It is more convenient to use the following equivalent 
representation 
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(") = P(R 2 ) (") + P(R 2 ) ( "' ) (22) 

\u>/ T \W f \u t -2uwj 

where P and Pi are polynomials with scalar constnat coefficients and Pq 

mentioned above may taken as Po(R) = P(R 2 ) R + Pi(R 2 )- 

Note that one could not apply immediately the Proposition 2.1' to this 



because the coefficient of dt in the representation (21) is not diagonal. On the 
other hand the operator R 2 has scalar leading part. First we will prove that 
if the symmetry (|22 ) admits an invariant surface then Pi in this equation 



vanishes. Let us take the invariant surface as u = q(w). Suppose that the 
function q(w) is differentiable at some point w = wq. Linearizing q around 
the point w (or as w — > w ) we obtain 

u - q(w ) = q'(w ) (w - w ) + o(w - w ). 

It follows from (|2l|) that in this case R 2 reduces to a scalar operator: R 2 — > 
(d t — w + q 2 (w )) I as w — > w , where I is the unit matrix. Thus in the 
linear approximation the Eq.(p2|) takes the form 



( U ) =P(d t -w + q\w )) ( U j +P 1 (d t -w + q 2 (w )) ( W j (23) 

where now P(d t —wo+q 2 (wo)) and P\(d t — wo+q 2 (wo)) are scalar operators. It 
is clear that the linearized equation is consistent with the linearized boundary 
condition u — q(w ) = q'(w ) (w — w ) , provided Pi = 0. Supposing that the 
equation ( p2|) is compatible with the constraint w = c where c is a constant 
and then linearizing about the point (u — 0, w — c) one can easily obtain 
that Pi vanishes in this case also. 
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It is evident now that in Proposition 2.1' one should put uq = 2, because 
R 2 = I dt + .... With this choice the constraint p(u, w) describes an invariant 
surface for the following system 

C),- R, C). (24) 

which is exactly the coupled Burgers type integrable system (see [8]) 

U T = U tt + 2(W + U 2 )u t 

(2b) 

w T = w tt + 2u 2 + 2(w + u 2 ) w t 

It is straightforward to show that the above system fl25|) is compatible with 
the constraint p(u, w) = only if p = w + u 2 + C\ u + c 2 or u = const. 

The above uniqeness proof of the boundary condition p = w+u 2 + Ci u+c 2 
can be more easily shown if we use a new property of the Burger's hierarchy. 
We have the following propositon: 

Propositon 2.3: The function u(t, x, r n ) (n > 1) satisfy infinitely many 
Burgers like equations 

u ,n, n ~ u tT2i+2 = -2 u^ D- 1 w Ti (26) 

Here i = —1,0,1,2,.... Burgers equation corresponds to i — — 1 , r_i = x 
, and To = t. All u n for % > —1 correspond to higher symmetries. It is 
straightforward to determine the even numbered symmetries of the Burgers 
equation from ( p6|) . It is very interesting that u satisfies the Burgers like 
equations with respect to the variables (r, , r 2 j + 2) for all i = — 1, 0, 1, 2, ... . 

The proof of this proposition depend crucially on definition of the higher 
symmetries of the Burgers equation. They are defined through the equation 
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u Tn = R n+1 u x (27) 

where R is the recursion operator given in eq.([TT]) and n > — 1. Eq.fl2"TD can 
also be written as u Tn = Ru Tn _ 1 . Differentiating this equation once by r n 
and using (|27|) one arrives at (f26|). 

If we let the most general boundary condition of the formp = f(u, u x ) = 
at x = Xq and take Tj and t 2 j +2 derivatives (for i > 0) of the function p and 
use the equation fl26|) we obtain 



Jm^ J,u,-u + Ju J,u x ,u x £ J t u x ^ J,u J,u x Ju,u x U \^°) 

Letting u = x\ and u x + u 2 + c\ u + c 2 = X2 then eq.(|28|) becomes 



J,x 2 J,xi,xi + J,xi J ,^ 2 ,x 2 ^ J ,x\ J,%2 J,xi,x 2 U, l^yj 

Assuming f X2 ^ and letting q = f, xi /f,x 2 we find that 

q,x 1 =qq,x 2 (30) 

This is a very simple equation and its general solution can be found. We 
shall not follow this direction to determine /(xi,x 2 ) rather change the form 
of equation p(u, u x ) = at x = x . This equation (in principle) implies either 

a) u x = h(u) 

which implies / = u x — h(u) at x = Xq . Or 

b) u = g(u x ) 

which implies / = u — g(u x ) at x = Xq It is now very easy to show that the 
cases a and b when the corresponding /'s are inserted in fl28|) we respectively 
obtain 
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a) h" + 2 = 

which implies u x + u 2 + C\ u + C2 = at x = x$ 

b) g" + 2 {g'f = 

which implies u = constant (for g = 0) and a special case of a (for g ^0). 
Hence we found all possible bounday conditions. 

Remark 2.3: On the invariant surface p(u, w) = the system ( p5|) turns 
into the Burgers like equation w T = Uu~ 2(ci M + C2) w* which is also integrable 
[5]. 

3 Application to Other Partial Differential 
Equations 

In this section we shall apply our method to obtain compatible boundary 
conditions of some nonlinear partial differential equations. Let us start with 
the following system of equations 

u t = u 2 + 2u 2 v 

—v t = V 2 + 2u V 
Letting v = u*, t — > it the above system becomes the well known nonlinear 
Schrodinger equation. Suppose that it admits a boundary condition of the 
following form 

u x \ x =o = p 1 (u,v) , v x \ x=0 = p 2 (u,v) (32) 

compatible with the fourth order symmetry. It means that the constraint 
( |32"D defines an invariant surface for this symmetry, presented as a system of 
four equations with two independent variables 
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u T = Uu — 2u 2 v t — 4:U v\ u\ + 2 v u\ — 2 u 3 v 2 , 

v T = —v tt — 2v 2 u t + Avu\ V\ — 2uv 2 + 2v 3 u 2 , 

ui,r = Ui t tt — 2u 2 vit — 2u\v\ — 6 u 2 v 2 u\ — 4 u v\ u t + 4 v u\ u t + 4 v v? vi, 

V\ jT = —Vijt — 2 v 2 ui tt + 2 v\ u\ + 6 v\ v 2 u 2 — 4 v u\ v t + 4 u v \ v t — 4 v 3 u u\ 

(33) 

One can check that the system ([33]) is compatible with the constraint U\ = 

p 1 (u,v),Vi = p 2 (u,v) only if p 1 = cu and p 2 = cv. Since the system (|33|) is 

of the form 

{u,Ux,v,v x )T = R 2 ^,^,^^)^ (34) 

hence it follows from the Proposition 2.1 that the constraints U\ — cu, V\ — 
cv are compatible with every symmetry of even order. So the boundary 
conditions u x \ x= q = cu , v x \ x= q = cv are compatible with such symmetries. 
Analytical properties of this boundary value problem are studied previously 
(see [2], [9-10]) by means of the inverse scattering method. 

Remark 3.1: On the invariant surface u\ = cu , v\ — cv the system 
( |33"D is reduced to a system of two equations: 

u T = u tt — 2u 2 v t — 2c 2 v 2 v — 2 u 3 v 2 , 

v T = —Vtt -2v 2 u t + 2c 2 v 2 u + 2 v 3 u 2 , 

The integrability of these equations is shown in [4] (see p. 175). Under a 
suitable change of variables in it this system of two equations becomes the 
famous derivative nonlinear Schrodinger equation. 

Among the nonlinear integrable equations the Harry-Dym equation 
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u t + u A u 3 = (35) 

is of special interest. It is not quasilinear and because of this reason its 
analytical properties are not typical. Using the symmetry approach we find 
a boundary condition of the form 

p(u,u x ,u 2 ) =0, (36) 

compatible with Harry-Dym equation. One has to notice that because of non- 
quasilinearity of ([35]) the transformation from the standard set of variables 
u, Mi, m 2 , u 3 ... to u, Mi, «2, u t , «i,t, M2,t, ••• is not regular. For instance u 3 = — ^. 
It has singular surface given by the equation u = 0. So one should examine 
this surface separately. Since the Harry-Dym equation ( |35D as well as its 
higher order symmetries possesses the reflection symmetry x — > —x, u — > 
—u,t — > t the trivial boundary condition w(£, 0) = is consistent with the 
integrability 

Suppose that the boundary value problem (|35|) and ( p6[) is compatible 



with the ninth order symmtery u T = u 9 Ug + .... It means that the constraint 
p(u, v, w) is consistent with following system of equations , equivalent to the 
ninth symmetry 

u T = /i 

v T = h (37) 

w T = h 
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where v = u x , w = u xx and (f u f 2 , f 3 ) T = R 3 (wt, v t , w t ) T , 

/ uw + u t d~j~ l w —uv — u t dt v u 2 + u t d^u \ 

R= - dt + v w — ^ + v t dt w —v 2 — v t d^v uv + vtdfu 

V w 2 + w t d^ l w -d t — vw — ^ — w t dt l v uw + w t d^ l u J 

The explicit expressions for f'2, fa are very long. Hence we give the explicit 
form only for the function f\: 

h = -uut + Zu tt ut\ - ~u tt u x h - §3- + §mti, tt /i+ 

(00) 
+\uu 1)t h t - ^uh 2 h t - ^h 3 u t - \uiu t h t . 

Where h = lu^u — u\. Here one has two choices for the rank of the 
equation fl36|). It is either one or two. The first choice does not lead to any 
regular invariant surface. The second gives 

1 , c 2 u . . 

u x \x=a = cu , u xx \ x= q = — — (39) 

Remark 3.3 : On the invariant surface v = cu, w = c 2 u/2 the first equation 
in the system (^) takes the form 

u T = -u m + 3u t u tt /u - 3u 3 u 2 /2 (40) 

equivalent to the MKdV equation. 

Since the symmetry under consideration is of the form u T = R 3 u x where 
R = u 3 D 3 u D^ 1 -5 the recursion operator for the Harry-Dym equation (see 
[11]) , the Propositon 2.1 implies the following 

Proposition 3.1: The boundary value problem ( |3"5| ) and fl3"6"| ) is compati- 
ble with every symmetry of the form u T = L(R 3 ) u x , where L is a polynomial 
with scalar constant coefficients. 
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The Korteweg de Vries equation u t = u xxx + 6u\u admits a recursion 
operator R = D 2 + 4w + 2u\D~ x which may be represented in the form: 

/ Au + Uvd^u l + 2vdt 1 \ 

R= d t + !2wdt l u -2u 2w d^ 1 

\2w+ 12(u t - 6uv) d^u d t - 2v -2u + 2{u t - 6uv) d^ 1 ) 

It is not difficult to show that the system of equations (u,v,w) T = H 3 (u,v,w) t 
admits an invariant surface u = 0, w = on which the equation turns into 
the MKdV equation. It means that the boundary condition u(t,x = 0) = 0, 
u xx (t, x = 0) = is compatible with all symmetries of the form u T = R 3n u x . 
Similarly, the MKdV equation u t = u xxx + Qu 2 u x is compatible with the 
boundary condition u(t, x — 0) = 0, u x (t, x — 0) = 0. 

4 Application to Discrete Chains 

Consider an integrable nonlinear chain of the form 

u t {n) = f(u{n- l),u(n),u(n+l)) (41) 

with unknown function u = u(n,t) depending on integer n and real t. The 
natural set of dynamical variables serving the hierarchy of higher symmetries 
for the chain is the set w(0), m(±1), m(±2), ... . However, it is more convenient 
for our aim to use the following unusual one, consisting of the variables 
u(0),u(l) and all their t-derivatives. Transformations of these sets to each 
other are given by the equation (|42|) itself and its differential consequences. 
In terms of new basic variables every higher order symmetry of this chain 

u T {n) = g{u{n — m),u(n — m — 1), ...u{n + m)), (42) 
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could be presented as a system of two partial differential equations 

v T = GiO, w, i>i, w ± ..., v s , w s ), 
w T = G 2 (v,w,vi,wi...,v s ,w s ), 

where v = u(0, t,r), w = u(l, t,r), v { = f~, Wi = f~. 
Prescribe some boundary condition of the form 

u(0)=p(u(l),u(2),...u(k)) (44) 

to the equation ( f4"TD to hold for all moments t. We shall call the boundary 
value problem (f41~D, (f4"4|) consistent with the symmetry (|4^ ) if the constraint 



4]) defines an invariant surface for the system (^3]). Note that interconnec 



tion between the hierarchies of the commuting discrete chains and integrable 
partial differential equations is well-known (see survey [4]). An illustrative 
example of such a kind connection is related to the famous Volterra chain 

u t {n) = u{n){u{n + 1) — u{n — 1)) (45) 

for which the next symmetry 

u T (n) = u{n)u{n + l){u{n) + u{n + 1) + u{n + 2)) — 
u(n)u(n — l)(u(n) + u{n — 1) + u{n — 2)) 
which might be represented as ([4], p. 123) 



Vt + % = (?VW + V 2 ) 4 , 

w T - w tt = (2vw + w 2 ) t 
under the substitution w(0) = v, w(l) = w, u(— 1) = w — a , 



(46) 



m(2) = v + — , u(— 2) = v ^r — -. Moreover, the full hierarchy of the 



Volterra chain is completely described by the hierarchy of the last system. 
According to the definition above the boundary value problem (f|4|), (fi5l) 
will be consistent with a symmetry of the Volterra chain if the constraint 



(f44|) describes an invariant surface for the same symmetry, represented as a 
system of partial differential equations. Let us examine invariant surfaces of 
the following system of partial differential equations 

v T = v ttt + (3vH 2 - 3v t H - 2v 3 ) t , 
w T = w ttt + (3wH 2 + 3w t H - 2w 3 ) t , 

H = v + w, which is exactly the higher order symmetry for the Volterra chain 
loT) of the form 



u T {n) = u{n)u{n + l){u{n + 2)u{n + 3) + u{n)u{n + 2) + u{n)u{n — 1) + 

u 2 {n) + 2u(n + l)u(n + 2) + u 2 (n + 2) + 2u{n)u{n + 1) + 

u 2 (n + 1)) — u(n)u(n — l)(u(n)u(n + 1) + u(n)u(n — 2) + u(n — 2)u(n — 3) + 

u 2 (n - 2) + 2u(n)u(n - 1) + u 2 (n) + 2u(n - l)u(n - 2) + u 2 (n - 1)). 

It is easy to check that the only invariant surface of the form v = const ad- 
missible by the system ( f44|) is v = 0. The corresponding boundary condition 
w(0) = is well studied (see [12], [13]). 

Remark 4.1: On the invariant surface v = the system (pE?D reduces to 
the scalar equation 

w T = w m + 3w tt w + 3w 2 + 3w t w 2 , 

which is nothing else but the next symmetry of the Burgers equation. More- 
over, the constraint is compatible with every generalized polynomial symme- 
try. On the invariant surface they are all reduced to the symmetries of the 
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Burgers equation. It is evident for instance, that the system fl4"T|) turns into 
the Burgers equation itself. 

Suppose now that v = p(w). Then one obtains that p(w) = —w. It gives 
rise to a boundary u(0) = — u(l) compatible with the Volterra chain (see 

[14])- 

Remark 4.2: Under the constraint v = —w the system (^) turns into 
the Modified KdV equation 

V T = Vm + Qv 2 V t . 

It is not difficult to show that there is no any invariant surface of the 
form v = p(w,w t ) such that J^- = admissible with the system fl37p. 

For the case Vt = p(v,w,Wt) calculations become very long so that here 
we utilized Matematica 2.1 (we thank G.Alekseev for his help with this cal- 
culations). Here p has a form p = —w t + 2v{v + w) which produces the 
boundary condition 

u(— 1) = —u(0) —u(l)—u(2). The slight difference with the (f44|) is overcomed 
by the simple shift of the discrete variable n. 

Using the proposition 2.1 it is easy to check that the invariant surface 
Vt = —Wt + 2v(v + w) is compatible with every odd order polynomial gener- 
alized symmetry of the system (|4"?D . It means that the boundary condition 
u{— 1) = —u(0)—u(l)—u(2) is compatible with the corresponding symmetries 
of the Volterra chain. 

The well-known boundary condition u 2 (0) = 1 for the modified Volterra 
chain 

u t (n) = (1 — u 2 {n)){u{n + 1) — u(n — 1)) 
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defines the invariant surface v 2 = 1 for the following systems of equations 

(48) 
w T -w tt = 2{{l-w')v) t 

and 



v T + v tt = 2((1 - v 2 )w) t , 
w T -w tt = 2((1 -w 2 )v) t 



(49) 



v T + v m = 2(u(l - v 2 )(3u> 2 - 1) - 3uura t ) t , 
w r + u>ttt = 2(w(l - w 2 )(3v 2 - 1) + 3vww t )t 
which are equivalent to the next symmetries of this chain: 

Ur {n) = (l - u 2 (n))(D_ - D+)(l - u 2 (n))(D_ - D + )u{n) 

and 

u T {n) = (1 - u 2 (n))(D_ - D+)(l - M 2 (n))[(-Z^ - D 2 )u(n) + 

p+ + D_)(n 2 (n)M(n + 1) + u 2 (n)u(n - 1) + 2w(n - l)«(n)«(n + 1)], 

here D + , D_ are the shift operators: D + u(n) = u(n + l), D_u(n) = u(n— 1); 
v = u(0), w = u(l) and other variables u(n) are expressed through v,w and 
their i-derivatives by means the chain and its differential consequences. 

Remark 4.3: On the invariant surface v 2 = 1 the systems (0), (|49| ) are 
reduced to the Burgers equation and its third order symmetry. 

5 Condition of weak compatibility 

It is easy to notice that any symmetry of the equation (1.1) rewritten in terms 
of the non-standard set of the dynamical variables turns into the equation 
containing m — 1 extra variables Ui,v>2, ...,w m -i • For instance, the fourth 
order symmetry of the Burgers equation 
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u T = -u 4 + 4u 3 u + 10u 2 ui + 6u 2 u 2 + 12-u^-u + 4uiu 3 
takes the following form 

U T = U U + 2(W + U 2 )u t , 

where w — u\. To extend it to the closed form it is enough to add one more 
equation obtained from the above equation by the differentiation with respect 
to x and replacing u 2 = Ut — 2uw. This is the general rule for integrable 
equations: one has to add m — 1 more equations (to have a closed system of 
equations), expressing variables Ui T , 1 < % < m — 1 through dynamical ones. 
But from the other hand side one may consider the single symmetry equation 
alone and suppose the extra variables are expressed interms of u and its lower 
derivatives. Let us pose the question , for which choice of such expressions 
the symmetry under consideration turns into an integrable equation? As 
an example let us consider the Burgers equation How should we choose the 
dependence w = w(u), such that the equation u T = u tt + 2(iu + u 2 )u t would 
be integrable? The only choice is w = —u 2 + c±u + c 2 (see [6]). We will 
call the boundary conditions Ui = Ui(u), x = (obtained this way) for the 
equation (1.1) as weakly compatible with the symmetry if these constraints 
are chosen to satisfy the requirement above: i.e. the equation for the nth 
symmetry written down in terms of the introduced variables turns into some 
integrable equation after replacing Ui = Ui(u), uu = Ut^£, .... So in the above 
case of the Burgers equation only the condition w(u) = —u 2 + c±u + c 2 is 
weakly compatible with the fourth order symmetry. As the remarks given 
above indicate, the compatibility of the condition with a symmetry implies 
the weak compatibility with it , but not vice versa. However, we conjecture 
that if the boundary condition is weakly compatible with at least three higher 
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symmetries then the corresponding initial boundary value problem will be 
solvable by a suitable generalization of the inverse scattering method. 

The following example for the Harry-Dym equation (|35| ) seems to be 
intriguing. Let us represent the fifth order symmetry 

U T5 = —-U 3 (2U5U 2 + IOU4U1U + IOU3U2U + 5u3ul) 

in the form u T5 = \ (hu) t , where h = 2u 2 u — u\. Represent also the next two 
symmetries in the similar form: 



u T7 = u tt ui - -u t uiuh + -u t [3(h + u\) 2 - Au\{h + u\) + uf] - uu m + -u 2t uh 

L O O 

and u T9 = fi (see above the first equation of the system (pTf)). It is evident 
that for arbitrary function F = F{u) the constraint h — 0, U\ — F(u) is 
weakly consistent with fifth and ninth symmetries, because the former takes 
the trivial form u T5 = and the latter turns into the integrable equation 
([40D. The seventh order symmetry becomes u^ = (Su t )t, where S = F — uF'. 
Thus, if for instance, S = a = const or S = , u ] 0) 2 one will have the equation 
u^ = (Su t ) t , to be integrable (see [5], p. 129). Supposing S(u) = a one can 

2 2 

easily find that u\ = cu + a, u 2 = ^ + ac + f-. It leads to the following 

2 

boundary condition u x = cu+a, u xx = |j, x = for the Harry-Dym equation, 
which coincides with (p9|) if a = 0. In the case S = , u } g \2 to find F one has 
to integrate the ordinary differential equation F(u) — uF'(u) = S. 
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